Deterministic preparation of Dicke states of donor nuclear spins in silicon by cooperative pumping 
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For donor nuclear spins in silicon, we show how to deterministically prepare various symmetric and asymmet- 
ric Dicke states which span a complete basis of the many-body Hilbert space. The state preparation is realized 
by cooperative pumping of nuclear spins by coupled donor electrons, and the required controls are in situ to the 
prototype Kane proposal for quantum computation. This scheme only requires a sub-gigahertz donor exchange 
coupling which can be readily achieved without atomically precise donor placement, hence it offers a practical 
way to prepare multipartite entanglement of spins in silicon with current technology. All desired Dicke states 
appear as the steady state under various pumping scenarios and therefore the preparation is robust and does not 
require accurate temporal controls. Numerical simulations with realistic parameters show that Dicke states of 
10 — 20 qubits can be prepared with high fidelity in presence of decoherence and unwanted dynamics. 
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I. INTRODUCTION 

Because of the ultra-long quantum coherence time, electron 
and nuclear spins of shallow donors in silicon are ideal can- 
didates as information carrier in quantum computation. For 
^^P donor in isotope purified ^^Si, a nuclear spin coherence 
time exceeding 1 s and an electron spin coherence time ex- 
ceeding 10 ms were measured^ ^. In his seminal paper^, Kane 
proposed an architecture for scalable quantum computation 
using nuclear spins of gated ^^P donors in silicon as qubits. 
With the superb compatibility with the existing silicon tech- 
nology, Kane's proposal and its variants have stimulated ex- 
tensive interests in donor systems in silicon^ ^, and remark- 
able technological progresses have since been made in vari- 
ous aspects. Accurate positioning of the ^^P donor was made 
possible by controlled single-ion implantation^. Local elec- 
trical tuning of the hyperfine interaction between donor elec- 
tron and nuclear spins was demonstrated^^ An architecture 
that integrates single ^^P donor with a silicon SET was re- 
cently developed which enables hi^h -sensitivity readout and 
control of the donor electron spin^^ In the meantime, chal- 
lenges still remain, e.g. on realizing efficient donor electron 
exchange coupling which is also needed for mediating pair- 
wise interaction of nuclear spins Due to the interference 
between the degenerate valleys of the electron, the exchange 
coupling strength vastly oscillates with sub-nanometer varia- 
tion in donor position (e.g. between 1-100 GHz when donor 
separation is 10 nm) The chance to have an exchange 
coupling of > 30 GHz expected by the initial Kane proposal 
becomes random unless with atomically precise dopant place- 
ments which is beyond the technological capability today and 
in the near future. 

In this paper, we introduce a scheme to deterministically 
prepare Dicke states of donor nuclear spins in silicon. Dicke 
states in general refer to the common eigenstates of and 
with eigenvalues J( J + 1) and M respectively, where J is the 
collective spin of an ensemble of N spin-/ particles. Dicke 
states with J taking the maximum value NI are symmetric 
under permutation operations. When M 7^ ±J, symmet- 
ric Dicke states are a class of genuine multipartite entangled 
states which have been widely pursued in atomic and optical 



systems as important resources for quantum information pro- 
cessing^^ These states have the remarkable properties that 
the entanglement is robust against qubit loss and projective 
measurements on the states lead to various entangled states 
of a lower qubit number. When J takes value other than 
NI, the Dicke states are no longer symmetric under all per- 
mutation operations and will be referred here as asymmetric. 
Asymmetric Dicke states are also resources of genuine mul- 
tipartite entanglement which are less studied for the lack of 
preparation schemes. The ideal resource for optimal quan- 
tum telecoloning algorithm is one such state with J = 0.1^ 
Preparation of asymmetric Dicke states also makes it possi- 
ble to access the decoherence free subsystems in the presence 
of collective decoherencJ^. The scheme we propose here 
can deterministically access all Dicke states (symmetric and 
asymmetric) that span a complete basis for the Hilbert space 
of N spins. This direct access to the Dicke-states basis can be 
an important complement to the circuit model quantum infor- 
mation processing, since the entanglement of these collective 
states can not be achieved in a simple way by pairwise inter- 
action 

Our scheme is based on cooperative pumping of nuclear 
spins by the coupled donor electrons, and the required con- 
trols are in situ to the prototype Kane proposal: (i) initializa- 
tion of electron spin to its ground state in magnetic field; (ii) 
ac electrical control of the hyperfine coupling of the donor (A- 
gate); (iii) on and off switching of exchange coupling between 
neighboring electrons (J-gate). Remarkably, our scheme only 
needs a sub-gigahertz exchange coupling, which can be sat- 
isfied for almost all donor pairs with separation ~ 10 nm. 
Hence it provides a practical way for generating the critical re- 
source of multipartite entanglement of spins in silicon, which 
can tolerate the exchange oscillation problem and be realized 
within the current technology. Our scheme is a significant ex- 
ample of the conceptually new approach of dissipative quan- 
tum state preparation with the advantages of robustness and 
no need for accurate temporal controls as compared to con- 
ventional state preparation by coherent evolution, as the de- 
sired Dicke states all appear as the unique steady state under 
the various pumping scenarios. Numerical simulation with 
realistic parameters for ^^P donors shows that Dicke states of 
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10 — 20 qubits can be prepared with high fideUty in the pres- 
ence of decoherence and unwanted dynamics. The scheme 
also appHes to other donor systems with larger nuclear spins 
such as ^^^Bi in silicoi#^^. 



II. CONTROL SCHEME 

Fig. 1(a) schematically illustrates the Kane architecture 
where shallow donors are embedded under patterned elec- 
trodes. The A-gate tunes the hyperfine interaction ^hf = 
^n^n ' with &n and In being respectively the elec- 
tron spin and nuclear spin of the nth donor. The coupling 
strength an is proportional to the electron density at the donor 
nucleus site and hence is a function of the voltage applied to 
the A-gate which pulls the electron wavefunction away from 
the nucleus. With each A-gate independently controlled with 
voltage Vn = Vn,o + ^Vn{t) cos{ujt + (j)n), the hyperfine in- 
teraction becomes: 

^hf = ^ CLn^oO-n ' In + ^^^^ COs{ujt + 0n)o-n ' In- 

n 

(1) 

Where the first term is the static hyperfine coupling where 
an,o/27r = a(Fn,o)/27r ^ 60 MHz, and the second term is 
the ac hyperfine coupling from the voltage modulation. In a 
strong magnetic field, the off-diagonal part of the static hyper- 
fine is far off resonance, hence it only results in small shifts 
of energy levels which are neglibile here. However, the off- 
diagonal part of the ac hyperfine can efficiently pump nuclear 
spin polarization when the modulation frequency is reso- 
nant to the electron-nuclear flip-flop transitioiP^M^. 

The exchange coupling between neighboring donor elec- 
trons is tuned by the J-gate which can separate the electrons. 
When the electrons are decoupled, nuclear spins of different 
donors can be independently pumped to the fully polarized 
state (i.e. nuclear spin initialization). When the electrons are 
coupled, nuclear spins can be cooperatively flipped by collec- 
tive raising/lowering operators. We have previously proposed 
using such pumping to probabilistic prepare singlets of nu- 
clear spind^^^. 

In general, for coupled donor electrons in uniform magnetic 
field along z direction, one can always find a non-degenerate 
one- spin-excitation eigenstate |eo) = Xln^^^nl^)' with 
\g) = \ ^ ■ ■ ■ \) denoting the ground state. With the ac control 
frequency uj tuned near resonance with the transition between 
1^) and |eo), we can neglect all other far detuned electron spin 
resonances. Dropping non- secular terms, the full Hamilto- 
nian of the electron nuclear spin system can then be written as 
H = Hq + He where: 

^0 = a;e|eo)(eo| +cjAr^/^ (2a) 

n 

- 15) (ffl aoin ^ |eo) (eo | aol^ , 

n n 

He = e*"*|eo)(<?| ^l^nWe*^"/- + hx. . (2b) 

n 
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FIG. 1: (a) Donor spins in silicon controlled by patterned electrodes. 
Pumping of nuclear spins is realized through ac voltage control of 
the A gates, (b) Energy level scheme for the exchange coupled elec- 
trons in magnetic field along z direction. The transition between the 
ground state \g) and a non-degenerate one-spin excitation state |eo) 
is used, (c) and (d) Pumping of nuclear spins by ac hyperfine cou- 
pling He (double head arrows), assisted by the electron spin decay 
channel from |eo) to |^) (wavy lines). The central frequency of He is 
tuned in resonance with the energy separation between |^) | J, M + 1) 
and |eo) I J^ M) which is a quantity dependent on M. (c) When the 
ac voltage control uses the first scenario for phases (see text), popula- 
tion is deterministically transferred from the collective nuclear state 
I J, M + 1) to I J, M). (d) When the ac voltage control uses the 
second scenario for phases (see text), population is transfered from 
I J, M + 1) to \J', M) where J' - J = or ±1. 



uje is the electron resonance frequency between |eo) and |^) in 
the magnetic field, ujn is the nuclear Zeeman frequency, and 
Vtn = §^^Vnan' The dc voltage of each A gate is set such 
that an,o = CLo- He induces electron-nuclear flip-flop where 
the nuclear spin flip is in a cooperative form determined by 
the phases and amplitudes SVn of the ac controls. For sim- 
plicity, we consider hereafter a uniformly coupled Heisenberg 
ring with the eigenstate |eo) = ^^{—l)^a^\g) which is 
gapped from other one- spin-excitation states by A (Fig. 1(b)). 
Then we have Qnit) = (-l)^l^(t) where n = ^§v^V. 

We consider two scenarios for the phases (j)n- In the first 
scenario, (— )^e*'^^ = 1 for all donors and hence the electron- 
nuclear flip-flop term becomes: He = l^(t)[e*^^|eo)(^| J~ + 
h.c.]. Here = J^±ijy where J = In is the collective 
spin of all nuclei. In the second scenario, (— )"^e*'^^ = 1 for a 
set of donors (referred as group A) while (— )^e*^^ = —1 for 
the rest (referred as group B), and hence the electron-nuclear 
flip-flop term is of the form: He = n{t)[e''^^\eo){g\Cj2 - 
Jb) + h.c], where Ja and Jb are the collective spin of nuclei 
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FIG. 2: (Color online) Preparation of the symmetric Dicke states of 
20 nuclear spin qubits initially on the fully polarized state. For the 
numerical simulation in (a-b), the amplitude Q of the ac control is 
held as a constant while the central frequency cj steps down after a 
finite interval. Symmetric Dicke states | J = 10, M) with different 
M are obtained sequentially with a probability > 99.8%. In (c-d), 
the control sequence Q{t) is optimized for speed, and each Dicke 
state is obtained with a probability > 99.5%. r/27r = 60 kHz. 



in group A and group B respectively, and Ja ^3b = J- Dicke 
states here refer to the common eigenstates of and with 
eigenvalues J{J ^ 1) and M respectively. Consider first the 
subspace with Ja = tiaI and Jb = {N — ua)! where ua and 
N — ua are the number of nuclei (of spin-/) in group A and 
B respectively, ja and jb are conserved quantum numbers in 
the dynamics. Dicke states in this subspace can be uniquely 
specified as | J, M). 

The first non-trivial element for preparing Dicke states 
is the deterministic population transfer from |J, M + 1) to 
I J, M). Only the first scenario for the phase control is needed 
where He conserves the quantum number J. We assume the 
population on electron excited state |eo) can be efficiently 
dumped to ground state with rate F, which can be realized, e.g. 
via the tunneling process between the donor and the SET is- 
land at low temperature^^ If the ac control frequency uj is in 
resonance with the transition |^) | J, M + 1) ^ | eo) | J, M) , the 
population on \g)\J^ M + 1) is pumped one-way to \g)\J-, M) 
via the Raman-type process [Fig. 1(c)]. For efficient initial- 
ization of the electron spin on the ground state, the electron 
Zeeman energy shall be large as compared to the tempera- 
ture. With a typical temperature of ~ 100 mK used in the 
experiments^, a magnetic field > 0.3 T is required, which 
corresponds to ac modulation of the A-gate voltage in the fre- 
quency range > 10 GHz. 

There exists two unwanted couplings by the electron- 
nuclear flip-flop: (1) the coupling to other electron one-spin- 
excitation states which are detuned at least by A [Fig. 1(b)]; 
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FIG. 3: (Color online) Preparation of asymmetric Dicke state | J = 
7, M = 7) in the subspace {ja — ^^Jb = 6) for 20 nuclear spin 
qubits initially on the fully polarized state, (a) Population flows by 
the cooperative pumping, (b) and (c) show numerical simulation of 
the state preparation, (b) Populations on the various Dicke states: 
I J = 10, M = 10) (dashed blue); | J = 9, M = 9) (dotted green); 
I J = 8, M = 8) (dotted red). The solid curve gives the popula- 
tion on the target state | J = 7, M = 7) (or its time reversal state 
in the shaded interval), and its final population is 99.5%. (c) The 
magnitude Q and central frequency uj for the ac tuning of hyperfine 
interaction. In the unshaded (shaded) interval, the first (second) sce- 
nario for the phases is applied (see text). The vertical arrows indicate 
timing of global tt flip of all nuclear spins. r/27r = 60 kHz. 



(2) the out-coupling of the final state \g)\J, M) to \eo)\J, M - 
1) which is also detuned by jj^ao. For the latter, we note that 
for any pair of states |^) | J, M+ 1) and | eo) | J^ M) coupled by 
He, the resonant frequency is — at + «o + ^<^o, namely, 
the pumping is M-selective with proper choice of the ac con- 
trol frequency. When the detuning S = min{A, j^ao} ^ T, 
these unwanted couplings cause a negligible leakage oc (^)^. 

With the nuclear spins initialized on the unentangled po- 
larized state I J = NI^ M = NI), the above simple control 
can deterministically prepare all symmetric Dicke states. Ex- 
amples for preparing symmetric Dicke states of N = 20^^P 
nuclear spins are simulated using the master equation 



p = i[p,H] - - (|eo)(eo|p + p|eo)(eo| 



2|^)(eo|p|eo)(^|). 

(3) 



The results are shown in Fig. 2. Clearly, the state preparation 
has the advantage that it is insensitive to the shape and area of 
the control pulse n{t). We take ao/27r = 60 MHz, r/27r = 
60 kHz, and assume the donor electron exchange of 0.5 GHz 
which leads to A/27r ~ 25 MHz. The preparation takes a time 
tp ~ 0(10) fis. In comparison, a single two-qubit nuclear 
spin gate mediated by such a small donor exchange would 
take ~ 1 ms^^. The probability to obtain each Dicke state 
is nearly unit (> 99.5%), and the imperfection is caused by 
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FIG. 4: Preparation of Dicke states \J, {ji^k}) which spans a 
complete basis for the nuclear spin qubits. (a) The quantum numbers 
ji^k are the collective spin of subsets of qubits. (b) Concatenated 
preparation process. In step /, Dicke states are parallelly prepared in 
each 2^-qubit subset with the collective spin being the specified value 
ji^k- The quantum numbers ji<iQ,k are all conserved in step Iq. Red 
solid (black dashed) lines indicate the donor electron exchange 'on' 
('off') when nuclear spins are pumped. 



the off -resonance couplings which maybe further reduced by 
using smaller Vt with the cost of longer preparation time. We 
further note that the reverse population transfer from | J, M) 
to I J, M + 1) can be realized if global tt flips of nuclear spins 
are applied before and after the above pumping process with 
the central frequency set at cj = cjg — ^at + fio — ^^^^^r^<^o- 
The second non-trivial element is to realize the asymmet- 
ric Dicke state |J = Jo,M = Jq) with general values of 
Jq. This requires the quantum number J to be changed which 
needs the control with the second scenario for ph ases. As dis- 
covered by the authors in our earlier wortf^*^, the inhomo- 
geneous collective operator — couples | J, M + 1) to 
the states | J + AJ, M) with the selection rule AJ = 0, ±1. 
Hence He can resonantly drive the one-way pumping from 
|J,M + 1) to |J + AJ,M) (see Fig. 1(d)). By repeating 
the above pumping for NI — Jq steps, the population can 
be transferred from the initial state \J = NI^M = NI) 
to the states |J > Jo,M = Jq) (see Fig. 3(a)). If the 
reverse pumping from |J, M = Jq) to |J, M = Jq + 1) 
(see proceeding paragraph) is also turned on, then all path- 
ways from the initial state will end up at the desired target 
state I J = Jo,M = Jq) where the population gets trapped 
(see Fig. 3(a)). This realizes the deterministic preparation of 
I J = Jo,M = Jo) with arbitrarily specified Jq. Fig. 3(b) 
shows the numerical simulation for preparing the Dicke state 
\J = 7,M = 7Ja = ^Jb = 6) of 20 ^^P nuclear spins, 
obtained with a probability of 99.5%. 

Combining these two elements, an arbitrary Dicke state 
can be deterministically prepared in the subspace defined by 
Ja = and Jb = {N — ua)!- Essentially, this is the real- 
ization of an arbitrary total spin eigenstate for the two collec- 
tive spins Ja and Jb • By concatenating this procedure, we can 
deterministically prepare any state in a complete Dicke states 
basis for N spins, denoted as | J, M, {ji,k}) where {ji^k} are 
the collective spin of subsets [Fig. 4(a)]. The concatenated 
procedures are illustrated in Fig. 4(b). 



TABLE I: Figure of merits for preparing Dicke states of 10 qubits in 
presence of decoherence. = | J = 5, M = 4) is the W-state. 
1^2) = I J = 5, M = 0) is the symmetric Dicke state with the most 
spin excitations, l^ps) = |J = 2,M = 2, j'a — §5 is — |) is an 
asymmetric Dicke state. F = {^plplip) is the fidelity of the density 
matrix with the target state at t = tp. We take r/27r = 50 kHz and 
is the peak value of Q{t) used. 
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III. EFFECTS OF POSSIBLE IMPERFECTIONS 

Here we analyze effects of various imperfections that may 
exist in realistic systems and give recipes on how to deal with 
these sources of errors. 

The first cause of error is the decoherence of the spin qubits. 
Dicke states are immune to decoherence in a collective form, 
but independent nuclear spin decoherence can cause a leakage 
out of the desired subspace. For a A/'-qubit state, the leakage 
is ~ tpNj where 7 is the nuclear spin decoherence rate and tp 
is the preparation timescale proportional to F"^. The longest 
nuclear spin coherence time reported for ^^P donor in sili- 
con is 1.75 second.*^ In Table. |l| we show the performance of 
preparing several exemplary Dicke states in presence of spin 
decoherence. Nuclear spin pure dephasing described by the 

Lindblad term E„ (iAp + plA " 2/^/^) is added 
to the master equation. The much slower nuclear spin relax- 
ation process is neglected^^. Because of the short preparation 
timescale, there is no visible effect from the decoherence if 
we take 7/27r = 0.1 Hz from the state-of-art measuremenffl 
Even with the much exaggerated decoherence rate 7/27r = 10 
Hz, the target states can still be obtained with high fidelity. 

In realistic systems, there also exist various defects in the 
sur rounding of the donors such as the interface P50 cen- 
^gj.g lio|37j[4orj]^g interplay between P50 centers and phosphorus 
electron spins is important in the presence of photo-excited 
electrons and holes. Such interplay can facilitates the elec- 
tric al det ection of paramagnetic resonance of the donor elec- 
^j.Qj{io|37Ji4ni jj^ absence of illumination as in our scenario, 
at temperature of RbT <C g/J^B, the P50 center as a mid- 
gap paramagnetic center can lead to statistical fluctuations of 
the electron-nuclear flip-flop resonances from donor to donor. 
Furthermore, there could be noises from the local electrostatic 
environment, e.g. possible charge traps in the gate oxidd^, 
which may affect the hyperfine coupling and the donor elec- 
tron exchange as well. If the charge hopping is faster than the 
state preparation timescale, the effect of this dynamics noise 
shall be similar to the effect of spin decoherence analyzed 
above. It is also possible to interlace the pumping control 
with sequence of tt pulses applied to the spin qubits for dy- 
namical decoupling from such noises. If charge hopping is 
much slower as compared to the state preparation timescale. 
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the noises are static which also resuU in the statistical fluctua- 
tion of the hyperfine coupling strength and exchange coupling 
strength. 

To deal with such statistical fluctuations of system param- 
eters, the system provides multi-fold individual tunability for 
each donor. The electron-nuclear flip-flop resonance depends 
on the static part of the hyperfine interaction (see Eq. (2a)) 
which can be independently controlled by each A-gate. Thus, 
in presence of inhomogeneous broadening of the electron- 
nuclear flip-flop resonances, the dc part of the A-gate volt- 
age can tune the static part of the hyperfine coupling for com- 
pensating the inhomogeneity after calibration. As the excited 
state for the electron-nuclear flip-flop has an intrinsic broad- 
ening of r, inhomogeneity in the flip-flop resonances is unim- 
portant if it is less than this intrinsic level broadening. 

The exchange couplings between donor electrons are in 
general inhomogeneous due to the exchange oscillation prob- 
lem and the possible effects of the aforementioned charge 
noises. This inhomogeneity manifests as nonuniform ampli- 
tudes an in the electron excited state |eo) = Xln^^^n 1^)- 
We note that the pumping operators that cooperatively flips 
the nuclear spin qubits are of the form: ^SVnOin^^^'^in 



(see Eq. pb])). While J-gate voltages can tune the exchange 
couplings and hence the values of a^, the ac modulation 
strength 5Vn and phase (j)n can further compensate the remain- 
ing inhomogeneity in a^. 

If only entanglement generation is of interest, our scheme 
can naturally cope with unknown systematic errors in the 
cooperative pumping. For controls aimed at the symmetric 
Dicke state (J")'^| J = N/2, M = N/2), the nuclear spins 
may be flipped instead by A~ = X]n(^ + ^nje*^""/" where 
On and r]n are respectively the unknown phase and amplitude 
errors. In presence of these errors, the steady state is still a 
definite pure state (A")^|J = 7V/2,M = with simi- 

lar multipartite entanglement. This feature makes the require- 
ment less stringent for implementing this scheme for entan- 
glement generation. 

Moreover, one advantage of our scheme is unwanted donor 



sites can be easily disconnected in the state preparation. There 
are two ways to do this. One is to turn off the coupling of 
such donor site with others by the J-gate control. Even if 
the exchange interaction couples all neighboring donor elec- 
tron spins into collective states, the set of nuclear spins being 
pumped can be further selected by the ac voltage controls. If 
the voltage applied to one of these donors has no ac compo- 
nent, its nuclear spin will not be pumped and is effectively 
decoupled from the rest in the preparation process. This ad- 
vantage allows one to select out the set of donors where the 
parameters are relatively uniform, and it also facilitates the 
realization of the concatenated preparation process described 
in Fig. 4. 

Finally, we give an estimate on the scale of the Dicke states 
that can be prepared using this approach. First, the detuning 
of the unwanted transitions is ~ -^ao which scales inversely 
with the qubit number N. To suppress the unwanted transi- 
tions, we require ^ F. Second, the leakage due to qubit 
decoherence is ~ tpNj, while the preparation time tp is de- 
termined by the smaller one of the modulation amplitude of 
the hyperfine interaction Q and the spin initialization rate F. 
In our simulation, we take F ~ 50 — 60 kHz and Q ~ 30 — 40 
kHz, both values in the range reported in experiment J ^ Q | ^ M 1 3 1 
Then N shall not exceed 100 by the first requirement men- 
tioned above. The upper bound for N set by the second re- 
quirement depends on the qubit decoherence rate 7. For given 
7, the fidelity for state preparation in a block of N qubits may 
be scaled from the numbers listed in Table I. For example, 
with the nuclear spin decoherence rate of 10 Hz (T2 ~ 15 
ms), preparation of asymmetric Dicke state {ips) for 20 qubits 
can have a fidelity over 0.96. Symmetric Dicke states may 
be prepared in a much larger scale, since they take much less 
steps by the cooperative pumping and tp is then considerably 
shorter. 

The work was supported by the Research Grant Coun- 
cil of Hong Kong under Grant No. HKU706711P and 
HKU8/CRF/11G. The authors acknowledge Andrea Morello 
and Wendi Li for helpful discussions. 



* wangyao@hkucc.hku.hk 

^ J. J. L. Morton, A. M. Tyryshkin, R. M. Brown, S. Shankar, B. W. 

Lovett, A. Ardavan, T. Schenkel, E. E. Haller, J. W. Ager, and S. 

A. Lyon, Nature 455, 1085 (2008). 
^ A. M. Tyryshkin, S. A. Lyon, A. V. Astashkin, and A. M. Raitsim- 

ring, Phys. Rev. B 68, 193207 (2003). 
^ E. Abe, A. M. Tyryshkin, S. Tojo, J. J. L. Morton, W. M. Witzel, 

A. Fujimoto, J. W. Ager, E. E. Haller, J. Isoya, S. A. Lyon, M. L. 
W. Thewalt, and K. M. Itoh, Phys. Rev. B 82, 121201(R) (2010). 

B. E. Kane, Nature 393, 133 (1998). 

^ R. Vrijen, E. Yablonovitch, K. Wang, H. W. Jiang, A. Balandin, 
V. Roychowdhury, T. Mor, and D. DiVincenzo, Phys. Rev. A 62, 
012306 (2000). 

^ A. J. Skinner, M. E. Davenport, and B. E. Kane, Phys. Rev. Lett. 

90, 087901 (2003). 
^ R. de Sousa, J. D. Delgado and S. Das Sarma, Phys. Rev. A 70, 

052304 (2004). 



^ J. J. L. Morton, arXiv: 0905. 4008 (2009). 

^ D. N. Jamieson, C. Yang, T. Hopf, S. M. Hearne, C. 1. Pakes, S. 
Prawer, M. Mitic, E. Gauja, S. E. Andresen, F. E. Hudson, A. S. 
Dzurak, and R. G. Clark, Appl. Phys. Lett. 86, 202101 (2005). 
L. Dreher, T. A. Hilker, A. Brandlmaier, S. T. B. Goennenwein, 
H. Huebl, M. Stutzmann, and M. S. Brandt, Phys. Rev. Lett. 106, 
037601 (2011). 

F. R. Bradbury, A. M. Tyryshkin, G. Sabouret, J. Bokor, T. 
Schenkel, and S. A. Lyon, Phys. Rev. Lett. 97, 176404 (2006). 
A. Morello, C. C. Escott, H. Huebl, L. H. Willems van Beveren, L. 
C. L. Hollenberg, D. N. Jamieson, A. S. Dzurak, and R. G. Clark, 
Phys. Rev. B 80, 081307 (2009). 

A. Morello, J. J. Pla, F A. Zwanenburg, K. W. Chan, K. Y. Tan, H. 
Huebl, M. Mottonen, C. D. Nugroho, C. Yang, J. A. van Donke- 
laar, A. D. C. Alves, D. N. Jamieson, C. C. Escott, L. C. L. Hol- 
lenberg, R. G Clark, and A. S. Dzurak, Nature 467, 687 (2010). 
S. Das Sarmaa, R. de Sousa, X. Hu, B. Koiller, Solid State Comm. 



6 



133, 737 (2005). 

R J. Ruess, L. Oberbeck, M. Y. Simmons, K. E. J. Goh, A. R. 
Hamilton, T. Hallam, S. R. Schofield, N. J. Curson, and R. G. 
Clark, Nano Lett. 4, 1969 (2004). 

H. Haffner, W. Hansel, C. F. Roos, J. Benhelm, D. Chek-al-kar, 
M. Chwalla, T. Korber, U. D. Rapol, M. Riebe, R O. Schmidt, C. 
Becher, O. Giiline, W. Dur, and R. Blatt, Nature 438, 643 (2005). 
R. Prevedel, G. Cronenberg, M. S. Tame, M. Paternostro, P. 
Walther, M. S. Kim, and A. Zeilinger, Phys. Rev. Lett. 103, 
020503 (2009). 

W. Wieczorek, R. Krischek, N. Kiesel, P. Michelberger, G. Toth, 
and H. Weinfurter, Phys. Rev. Lett. 103, 020504 (2009). 
M. Eibl, N. Kiesel, M. Bourennane, C. Kurtsiefer, H. Weinfurter, 
Phys. Rev. Lett. 92, 077901 (2004). 

L.-M. Duan and H. J. Kimble, Phys. Rev. Lett. 90, 253601 (2003). 
A. S. S0rensen and K. M0lmer, Phys. Rev. Lett. 91, 097905 
(2003). 

J. K. Stockton, R. van Handel, and H. Mabuchi, Phys. Rev. A 70, 
022106 (2004). 

C. Thiel, J. von Zanthier, T. Bastin, E. Solano, and G. S. Agarwal, 

Phys. Rev. Lett. 99, 193602 (2007). 

L. M. Duan, Phys. Rev. Lett. 107, 180502 (2011). 

M. Murao, D. Jonathan, M. B. Plenio, and V. Vedral, Phys. Rev. 

A 59, 156, (1999). 

J. Kempe, D. Bacon, D. A. Lidar, and K. B. Whaley, Phys. Rev. A 
63, 042307 (2001). 

R. E. George, W. Witzel, H. Riemann, N. V. Abrosimov, N. 



Notzel, M. L. W. Thewalt, and J. J. L. Morton, Phys. Rev. Lett. 
105, 067601 (2010). 

M. H. Mohammady, G. W. Morley, and T. S. Monteiro, Phys. Rev. 
Lett. 105, 067602 (2010). 

G. W. Morley, M. Warner, A. M. Stoneham, R T. Greenland, J. 
van Tol, C. W. M. Kay, and G. AeppH, Nat. Mat. 9, 725 (2010). 
M. S. Rudner and L. S. Levitov, Phys. Rev. Lett. 99, 246602 
(2007). 

E. A. Laird, C. Barthel, E. I. Rashba, C. M. Marcus, M. P. Hanson, 
and A. C. Gossard, Phys. Rev. Lett. 99, 246601 (2007). 
E. L Rashba, Phys. Rev. B 78, 195302 (2008). 
W. Yao, Phys. Rev. B 83, 201308 (2011). 

H. Y. Yu, Y. Luo, and W. Yao, Phys. Rev. A 84, 032337 (2011). 

C. D. Hill and H.-S. Goan, Phys. Rev. A 68, 012321 (2003). 

D. R. McCamey, J. V. Tol, G. W. Morley, and C. Boehme, Science 
330, 1652 (2010). 

A. R. Stegner, C. Boehme, H. Huebl, M. Stutzmann, K. Lips, and 
M. S. Brandt, Nat. Phys. 2, 835 (2006). 

H. Huebl, F. Hoehne, B. Grolik, A. R. Stegner, M. Stutzmann, and 
M. S. Brandt, Phys. Rev. Lett. 100, 177602 (2008). 

G. W. Morley, D. R. McCamey, H. A. Seipel, L.-C. Brunei, J. van 
Tol, and C. Boehme, Phys. Rev. Lett. 101, 207602 (2008). 

H. Morishita, L. S. Vlasenko, H. Tanaka, K. Semba, K. Sawano, 
Y. Shiraki, M. Eto, and K. M. Itoh, Phys. Rev. B 80, 205206 
(2009). 

Z. R. Gong and W. Yao, in preparation. 



